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SYMBOLS 


2a,  26, 2c 
u 

2^]  6, 
l^m^n 

^mnj  BnuDi 

m 

/(s^,2) 


Cartesian  co-ordinate  axes 
Length,  depth,  breadth  of  the  prism 
Poisson’s  ratio 

Normal  stress  components  in  the  x,y,z  directions 
Lengths  of  loaded  portion  in  the  y  and  z  directions 
Indicies  of  the  Fourier  series 
Fourier  coefficients 

Stress  distribution  function  on  the  faces  x  =  ±a 


1. 


INTRODUCTION 


The  stress  distribution  in  rectangular  prisms  subjected  to  certain  normal  end  forces 
on  its  two  opposite  faces  is  a  fundamental  problem  in  Applied  Mechanics.  The  stress 
distribution  near  the  loaded  surface  may  be  very  complex,  and  the  distribution  given  by 
elementary  theory  is  often  inaccurate  [1].  In  this  paper,  it  is  assumed  that  the  bulk  stress 
((Tz  -f  +  (Tx)  on  a  lateral  surface  of  a  rectangular  prism,  which  is  being  subjected  to 
an  arbitrary  losiding,  is  known  (i.e.  measurable).  Such  a  measurement  could  be  obtained 
using  advanced  thermal  emission  techniques.  The  aim  of  this  paper  is  to  reveal  how  the 
measured  values  of  bulk  stresses  csm  be  used  to  determine  the  actual  end  loads. 

A  method  for  calctdating  the  bulk  stresses  produced  by  a  self-equibbrating  end  load 
is  given  in  [1].  This  will  be  referred  to  as  the  Direct  solution.  In  the  Inverse  method  we 
will  determine  the  end  loading  from  a  known  distributior  of  bulk  stress  values  on  a  lateral 
surface.  In  this  case,  biJk  stresses  on  the  lateral  surface  were  simulated  using  the  Direct 
solution.  In  this  preliminary  “proof  of  concept”  work,  the  effect  of  experimental  noise  on 
the  thermal  signal  was  not  considered. 

This  technique  of  determining  the  actual  end  load  based  on  the  known  brdk  stress  field 
has  many  potential  applications,  i.e.  in  fatigue  tests  involving  rectangular  cross  sections 
to  determine  the  stress  distribution  inside  the  specimens.  Other  possible  applications 
include  the  determination  of  the  longitudinal  stress  distribution  around  a  planar  crack  at 
right  angles  to  the  sides  of  the  '  rism. 


2.  DIRECT  SOLUTION 

The  solution,  given  in  [1],  is  first  used  to  generate  the  bulk  stresses  on  a  lateral  surface 
for  the  case  of  a  rectangular  prism  subjected  to  end  forces  symmetric  in  x,y  and  z.  For 
the  rectangular  prism  shown  in  Figure  1,  in  the  absence  of  the  body  forces,  the  general 
solution  is  as  follows. 
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First  define  the  Galeikin  vector  F. 
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The  stress  components  are  then  related  to  F  by  the  following  equations. 

«r,  =  2(1  -  +  (•'V'  -  div  F,  (2) 

r,..(l-.)(|VF,  +  |Vf.)-^di,r.  (3) 


2 


> 


This  gives: 


=EE 


■AmnOf 

cosh  am-a 


[(1  +  Omn^  coth  am„a)  cosh  a^nX  -  a^nX  sinh  amnx] 


xrvKy  nirz  Bnio 

xcos— cos— +  EE^r:3 — 


T"  I  cosh^„(6 


2un^ir^  coshjSniK  4- 


X  [( 1  -  coth  /3n|h)  cosh  +  /3„,y  sinh  /9„iy]  I  cos  —  cos  — 


+ 


EE^ 


DlmC 


I  roth-llmC 

X  cosh  7jmz  +  7/m2  anh  Tlm^j 


—  j  [(1  -  7imCCOth7i„ 


lifx  miry 

cos - cos  — ; — 

a  b 


(4) 


_  _  Arnnamnir^ 

~  ^  ai,.iccoshamn<i 


m  n 

4"  (^mnX  sinh  CXfnnX 


[1  -  (2l>  +  OmnO  coth  am„a)]  cosh  Qm„I 
B„innb 


.  nirz  .  miry 

s>“ - an  — ; - /  >  „  ■  ^  . 

c  fc  n  I 


/3„(y  cosh/?„jy 


-  3nib  coth  0„ib  sinh  0niy 


nirz 


lirx 


Dimmirc 


lirx 

sin - cos - 1-  >  >  ■ 

c  a  ^ "  7)mhcosh7i„c 


7Jm2  cosh7(mZ  -  7lmCcoth7j„csinh7(„z 


I 

.  miry  lirx 

sin  — ; —  cos - 

b  a 


(5) 


Expressions  for  the  other  stress  components  can  be  written  in  a  similar  manner. 
With  reference  to  Figure  1,  the  following  boundary  conditions  apply. 
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All  the  shear  stress  boundary  conditions  are  exactly  satisfied  by  (5)  and  similar  expres¬ 
sions  for  and  The  normal  stress  boundary  conditions  are  approximately  satisfied 
by  equating  low  order  Fourier  cosine  terms  in  equation  (4),  and  in  similar  expressions  for 
cTy  and  (Tz.  Taking  the  double  Fourier  cosine  transforms  of  these  three  equations  (with 
their  approximated  normal  stress  boundary  conditions)  gives  the  expressions  (7), (8)  and 
(9). 
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where 
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where 
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Equations  (7), (8)  and  (9)  may  then  be  solved  simultaneously  for  the  unknown  Fourier 
coefficients.  The  stress  at  any  point  in  the  prism  may  then  be  calculated  by  substituting 
the  Fourier  coefficients  into  equations  (4), (5)  and  similar. 

Since  the  given  equations  are  valid  only  for  self  equilibrating  end  loads,  it  is  necessary 
to  modify  the  required  end  loading  f(y,z)  so  that  the  net  force  acting  on  each  end  face 
is  zero.  This  produces  a  modified  loading  distribution  fj[y,z)-  For  symmetric  loading, 


h(y,z)=  ny<^)--^ 


where  P  is  the  total  loswl  due  to  f(y,z),  i-e. 

r+4  r+c 


J  ‘  f(y,^)dydz 


(10) 


(11) 


The  previous  equations  calculate  the  stresses  due  to  the  modified  loading  /i(y,  z).  The 
stress  component  <2^,  calculated  &om  equation  (4),  is  a  Fourier  series  approximation  to 
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3.  INVERSE  SOLUTION 

The  main  aim  of  this  section  is  to  provide  an  Inverse  solution  to  the  problem,  which 
involves  calculation  of  the  end  stress  from  the  known  bulk  stress  distribution  {a^  + 
<Ty  +  <Ti)  on  one  lateral  surface  of  the  prism.  Consider  the  case  when  we  know  the  bulk 
stress  distribution  on  the  surface  y  ~  +b.  On  this  surface,  o,  =  0,  and  so  the  bulk  stress 
is  o*  +  Oj.  From  equation  (4)  we  obtain  the  following  expression  for  the  bulk  stress. 


{(Ti  +  <Tz 


'v=(i 


=  EE 


cosh  a„„a 


-  a™„isinha„„i  cos 


(1  +  QmnOCOth  OmnO)  COsh  Qmtir 

Bn,b 


TtlTZ  ^ 

+  EE 


n  / 


0liC^  cosh 


(Ittc  \  * 

- j  [( 1  —  0nlbcoih  0„ib)  cosh  d„ib 

DlmCi-ir 


+  0„lbsmh0„ib\ 


Tiirz  Ittx 

cos - cos -  +  >  > 

^  ^  l  m 

2 


yf„b^  roshT(„r 

2i/m^jr^  cosh  7j„2  +  1(1  -  7lmCcoth7imc) cosh  7jm2 

iirx  ~  ~  " 


+  7im^sinh7i„jzj 


cosh  'tlmC 


C08h7(„2  -  7(„rsinh7imZ  cos 


Ittx 


EE 


(1  +  7,„ccoth7/„c) 

-4,„na(-ir 


2vrr0ir^  cosh  a„,„7  + 


4“  Onin®  sinh  Onirt^] 


„  „  “mn^  “mnO 

(=?)■ 

TITT  Z  yt  ^ 

c  0^,0?  cosh  0.nlb 


1  -  Q„„ocothu,„na)co8ha„,„j- 
_ Bnih 

n  I 


2v0ir^  cosh  0nib  +  ( - )  [( 1  —  0nlb  coth  0nib)  cosh  /3n(6 


+  /3„i6sinh/3„(6) 


nx2  lirx 

cos - cos - 

c  a 


(12) 


6 


r 


> 


I 


1 


By  taking  the  Fourier  transform  of  this  equation,  as  shown  in  Appendix  A,  the  following 
results  are  obtained. 
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Equation  (13)  may  then  be  applied  for  each  value  of  n  and  I  to  give  a  set  of  linear 
simultaneous  equations  with  the  Fourier  coefficients  A,„„,  B„i  and  Di^  as  the  unknowns. 
Since  the  boundary  conditions  of  equation  (6)  still  apply,  equations  (7)  and  (8)  are  also 
used  in  conjunction  with  (13).  In  effect,  equation  (13)  replaces  equation  (9)  for  the  Inverse 
solution.  This  imposes  the  restriction  that  I  =  m,  so  that  the  number  of  linear  equations 
produced  by  (7)  (8)  and  (13)  equals  the  number  of  unknown  Fourier  coefficients.  The 
resulting  set  of  equations  is  then  solved  to  produce  a  set  of  Fourier  coefficients  Amm  B„i 
and  Dim  Inverse  problem.  Equation  (4)  is  then  used  to  obtain  the  applied  stress 
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4. 


RESULTS  AND  DISCUSSION 


4.1  General 

The  Direct  and  Inverse  programs  were  implemented  using  FORTRAN  77  on  an  Apollo 
DNIOOOO  computer.  A  variety  of  numerical  procedures  were  considered  for  performing 
the  integrations  required  by  equations  (9)  and  (13).  Initial  versions  of  the  programs  used 
either  trapezoidal  or  Simpson’s  rule.  It  was  found  that  these  methods  were  not  sufficiently 
accurate,  and  the  Gauss-Quadrature  method  was  finally  selected.  Reference  [2]  describes 
the  Gauss-Quadrature  technique  in  detail. 

Various  methods  were  also  used  for  the  solution  of  the  resulting  set  of  linear  simul¬ 
taneous  equations.  Gaussisin  elimination  with  full  pivoting  was  employed  in  the  final 
program. 

Both  the  Direct  and  Inverse  programs  were  written  with  the  aim  of  minimising  nu¬ 
merical  errors.  Double  precision  was  used  throughout  the  programs,  and  many  mathe¬ 
matical  functions  were  re-written  in  more  suitable  forms  to  reduce  errors  for  large  values 
of  1,  m  and  n. 


4.2  Direct  Results 

The  equations  given  in  the  previous  section  are  applicable  to  a  prism  with  an  arbitrary 
rectangular  cross  section.  However,  for  iUustrative  purposes,  a  square  cross  section  with 
a  square  loading  area  was  considered.  In  this  case,  b  ~  c  ~  2  and  k\  =  ki  =  0.5. 
With  reference  to  Figure  1,  a  stress  of  f(y,z)  =  16  was  applied  to  the  l-adcd  area,  with 
f(y,z)  =  0  outside  the  loaded  area.  Equations  (10)  and  (11)  must  first  be  used  to  produce 
a  self-equilibrating  loading  /i(v,z)  on  the  end  face. 
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From  (11)  the  total  load  is; 


P  =  y  f{y,z)dydz 
=  16(2ifei6)(2jfe2c)  +  0 
=  64/e|  1:26c 


Applying  (10)  gives  the  modified  loading  function. 


Hence, 


/i(y.2)  = 


=  /(k.-^)  - 


46c  J 

641;]  62  6c 


46c 


=  /(y,z)-4 


fl(y,  2)  =  16  —  4  =  +12  over  the  loading  area 
=  0  —  4  =  — 4  outside  the  loading  area 

This  applied  stress  field  is  shown  in  Figure  2. 

Sample  results  from  the  Direct  method  are  shown  in  Figure  3.  Plots  of  the  bulk  stress 
distribution  on  the  y  =  +6  surface  tind  the  end  stress  (r*  are  presented.  These  results 
were  obtsdned  using  seven  terms  of  the  Fourier  series  (i.e.  l,m,n  —  0,1,  ...7),  and  show 
good  convergence.  Shown  in  Figure  3a  is  the  surface  plot  of  (t*,  with  its  corresponding 
values  indicated  in  the  contour  plot,  Figure  3b.  These  plots  show  the  stress  values  ranging 
from  —5.0  to  +14.0  with  the  majority  between  —4.0  to  +12.0.  This  solution  is  in  good 
agreement  with  the  applied  stress  field  /i(y,  2).  However,  a  small  discrepancy  is  noticeable 
around  the  edges  of  the  surface.  The  bulk  stress  distribution  on  the  y  =  +6  surface  is 
shown  in  Figures  3c  smd  3d.  As  expected,  symmetries  occur  about  the  z  =  0  and  2  =  0 
axes  with  zero  stress  values  occuring  sJong  the  2  =  0  axis. 
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4.3 


Inverse  Results 


In  the  Inverse  problem,  the  bulk  stress  field,  generated  by  the  Direct  method,  was  used 
to  calculate  the  applied  end  stress.  For  comparison,  the  results  of  the  stress  distributions 
obtained  from  both  the  Direct  and  Inverse  methods,  at  different  Fourier  indicies,  are 
shown  in  Figure  6. 

Initial  program  versions  produced  good  Inverse  results  for  all  Fourier  indices  up  to 
the  fourth  term.  However,  at  higher  terms,  the  Inverse  solution  was  poor.  The  problem 
occurred  evaluating  the  double  integral  of  equation  (13),  using  Simpson’s  rule  with  up  to 
161  grid  points  in  each  direction.  A  Gauss-Quadrature  method  was  then  implemented,  in 
a  general  form,  so  that  the  number  of  Gauss  points  could  be  increased  easily  if  required. 
By  using  81  grid  points  in  the  x  and  z  directions,  good  results  were  obtained  for  up  to 
eight  Fourier  indicies  (i.e.  l,Tn,n  =  0,1,. ..8).  The  number  of  Gauss  points  used  between 
each  grid  point  is  indicated  in  Figure  4. 

It  is  evident  that  more  Gauss  points  must  be  used  when  solving  for  a  higher  number 
of  Fourier  coefficients.  Figure  5  shows  the  effect  of  increasing  the  number  of  Fourier 
coefficients  without  changing  the  number  of  Gauss  points  used  by  the  integration  routine. 
Figures  5a  and  5b  show  the  surface  and  contour  plots  produced  by  the  Direct  program  for 
l,m,n  =  0, 1,  ...9.  Figures  5c  and  5d  show  similar  plots  produced  by  the  Inverse  program. 
Clearly  the  agreement  between  the  two  solutions  is  not  as  good  as  for  the  plots  presented 
in  Figure  4.  Further  increases  in  the  number  of  Gauss  points  and  Fourier  coefficients 
were  not  warranted  for  preliminary  work  due  to  the  large  increase  in  computer  resources 
required. 

One  disadvantage  of  the  Gauss-Quadrature  method  is  that  the  bulk  stress  distribution 
is  required  at  points  which  are  not  evenly  spaced  over  the  surface  of  interest.  Since  a 
thermoelastic  (i.e.  SPATE)  measurement  is  likely  to  contain  stress  information  for  a 
number  of  evenly  spaced  points,  the  stress  values  at  the  location  of  the  Gauss  points 
would  first  need  to  be  interpolated  from  the  measured  values. 
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5. 


CONCLUSION 


This  work  has  shown  that  it  is  possible  to  predict  the  stress  distribution  in  a  rect¬ 
angular  prism  using  bulk  stress  measurements  &om  one  surface  only.  Analysis  has  been 
performed  for  the  case  where  the  bulk  stress  distribution  used  in  the  Inverse  method  was 
produced  by  the  Direct  method.  The  results  obtained  were  good  for  up  to  8  Fourier 
indices. 

Although  the  concept  has  been  successfully  proven,  a  considerable  amount  of  work 
will  be  required  to  extend  this  methodology  to  a  procedure  which  can  be  routinely  used 
for  the  analysis  of  thermoelastic  measurements  from  real  structures. 
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Fig.  3b 


Fig.  3d 


FIGURE  3.  DIRECT  SOLUTION  FOR  l,m,n  =  0,1,  ...7 


! 


Direct  Solution 


Inverse  Solution 


Fig  4n.  l,m,n  =  0, 1.  Two  point  integration. 


Direct  Solution  Inverse  Solution 

Fig  4b.  I,m,n  =  0, 1,2.  Two  point  integration. 


Direct  Solution 


Inverse  Sointioa 


Fig  4c.  i,m,n  =  0, Two  point  integration. 


> 


i 


Direct  Solution 


Inverse  Solution 


Fig  4g.  I,m,n  =  0,1. ..7.  Five  point  inlegiation. 


Fig  4h.  /,m,n  =  0, 1...8.  Six  point  integration. 


FIGURE  4.  DIRECT  AND  INVERSE  SOLUTIONS  FOR  VARYING  l,m,n 


r 


i. 


Fig.  6a.  Direct  Solution 


Pig.  6c.  Inverse  Solntion 


Fig.  5b.  Duect  Solution  Fig.  5d.  Inveiue  Solution 


FIGURE  5.  DIRECT  AND  INVERSE  SOLUTION  FOR  l,m,n  =  0,1,  ...9 


APPENDIX  A 


DERIVATION  OF  EQUATION  (13) 


From  Equation  (4) 


+  <72  =  Y'  y'  — — ((1  +  am„ocothQ„„o)cosha„„i  -  o„„ *  sinh  q„„ r] 
cosha™„a 


miry  nTi  r— ^  B„|6  2  2  ,  „  ( \  ^ 

X  cos  cos - +  2^  V  ^52— 2 - r3~i:  cosh^niy  +  — 

b  c  ■‘^■^B„,c’cosh;3„)6  [  \  a  J 

[(1  -  0nlb cothflnib)  cosh 0„iy  +  0„iysm}t0„iy]  cos  cos  — 

c  a 


+  W _ _ _  : 


cosh  7jfn-  + 


-)i„ccoth ->i„c) 


.  -21 

X  cosh-rimC  +  7imCsmh7/„zl  cos - cos—— 

a  0 


y^  'T'  +  7,„ccoth7(„c)cosh7i„z  -  7j„z  sinh  7/„z] 

cosh7i„c 

I  m 

Inx  mxy  ^V*'  ^mna  2  l  f  nirb\‘ 

X  cos - cos  — 7— +  >  >  ,  •  .. - r -  2i/m"n-coshatn„x  +  - 

a  62cosho,nna  \  c  / 


_ Atnn^ _ 

^tnn^ 


f/  t  V  I  •  1  ^  »Ti»ry  nnz 

X  1(1  -  a„„acotha„„o)coshQ«„x  +  q„„x sinh Qto„xJ  cos  — j—  cos - 

0  c 


v'  Y'  Bolt 

„  I  coshi3„it 


,  .  _  -  .  2.  1  '•’'■2  lirj 

X  coshBniy  +  BnlJfSinhBolJl)  cos  -  cos  - 

c  a 


2W^ir^coshl3„iy  +  (-^)  1(1  -  0nibcoHi0„ib) 


On  a  lateral  surface  of  y  =  6,  Oj,  =  0  and  so  the  bulk  stress  on  this  surface  becomes 
<7'fc  =  <7i  +  <72.  Putting  y  —  b  into  the  above  equation  yields  the  following  result. 
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"■i.  =  ^ - ^1(1  +  Omnacotha„„a)co5ha„„x  -  a„„xsinha„„Tjcos 

‘—‘^-1  cosha„„a  c 


+Ei: 


B„,6 


0l,c^  cosh  A,i5 


n  J 

X  cosh/3„|4  -!-  Ai|isinh/3„ii] 


/  Irr  \  ^ 

2*xn’«'’cosh/3„i6+  [(1  - /3„i6coth/3nil>) 


nrz  lirx 

cos - cos - 

c  a 


ZJlm< 

^  rL6’  < 


X  cosh  Tim  J  +  -Vim  X  sinh  ->)„  2 

DlmC(-l)’ 


2t/m  )r  cosh  Tim  z  + 
Ivx 


(t)  k 


1  -  'yimCCOth  -yimC) 


EE 


coshTimC 


1(1  +  TlmCCOth  TlmC)coshTlm2  -  Tln.2  sinh  Tim  2)  COS 


Ixx 


A„«a(-ir 


cosh  Ckmn® 


CO>)xa„^X  +  1(1  -  OmnOfOth  OmnO) 


X  COsham>.2  +  Omni  sinh  OmnX 

flnifc 


X  cosh/3„i6  +  /3„|6sinh/3„|61 


IT*  cosh i5„;6  +  [(1  _  t)„,6coth 


l^nrh) 


nit:  litx 

cos - cos  — 


(I) 

(II) 

(Mil 

(IV) 

(V) 

(VI) 


To  determine  the  Fourier  coefficients  A^ni  B^i  and  Dim  >•  is  necessary  to  take  the  Fourier 
transform  of  the  above  equation. 


f*" I  1  j  j 

/  /  (<r, +  <r,  )cos - cos - dxdx 

J-c  J-a  o  c 


-i:n 


(i)4-(n)  +  (iii)-»-(iv)-»-{V)-^(vi) 


pwi  qxz 
cos - cos - dx  dz 


Note  that  p  and  q  are  now  fixed  positive  integers. 


appendix  a 


when  n  =  0 
when  n  ^  0 


where 
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Taking  term  (II) 


/  /  (Ii)cos - cos - axdz 

J-c  J-.  “ 

B„b 


''  ^  ?  01,  c 


cosh  (3nlh 


2vtirir  coshf^ni^  4- 


(t)’« 


1  -  0„,bcolh0„,b) 


■  .  j  9’’'-  j 

cosh/3„|fc  + /3„i6sinh/3„|6]  /  cos - cos - dx  cos - cos - d: 

\J-a  “  “  J-c  c  c 

B„iab(l  +  S„i)(l  +  6^)  f2j,„2y2cosl,^_^,5  ^  — ^  ((1  -  l3„ib  colb  0„,b)  cosb  3„,b  +  /3„|6sinh)3„/fc] 

2i/n*ir’ +  1(1  - /3„l6coth/3„|fc)  + /3„|6tanh/3„i6] 


0l,c  cosh  0„ib 

Bnia6(l  +  ^oj){l  +  bon) 


01, C 


where 


iol  |q 


when  1  =  0 
when  1^0 

when  n  =  0 
when  n  ^  0 


Taking  term  (III) 


I"  cos  ?^cos^-^dxdz 


__  ^  AmC(  — 1)"* 

~  f  ^  cosh^tmC 


(t) 


(1  -  7lmCCOth  7l„c)  , 


r’  u  1^- X  ■  K 

j  cosh  7i„r  cos - a;  +  — j  ■yim  J  zstnn 

I'KX  P'XX 


qirz 

■yirr,  z  cos - az 


■r 


COS -  COS 

a 


-  dx 


Z  T  [{  +  (  ^)  ■  (1  -  7,™ccoth  7,„c)| 

^  7,’„6*cosh7imC  \  o  /  J 


27ltn(-l)’  SinhTlmC 

liL  +  (?»/<^)’ 


(? 


27in.c(-l)’  cosh  71,. 

iL  +  (9’f/f)^ 
(?*-/c)'^ 


V"  Am2oc(-l )"’*’(!  +#,/)  f _ J  J _ 

=  >  - rm - ; — ,  »  tanh  7imf  + 

^  7i™»’(Ti'„  +  (Wc)M  . 


I  —  1  (1  -  7/«iCCOlh7|mC) 


xtanh7,™c+(-)  -^-^^^tanh 


7lmC 
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Put  n  =  q,  re-arrange  and  simplify  to  give; 


=  1: 


4g|„  ac(  - 1  1  +  )  tanh  ->i„c 

iyf„  +  (nir/c)’) 


{lirb/a)^(nir/c)~ 

■yL  + 


where 


r  1  when  1  =  0 
I  0  when  I  /  0 


The  integrals  of  parts  (IV)  (V)  and  (VI)  were  derived  using  the  same  approach  as  shown 
in  (I)  (II)  and  (III). 


■  COS - dx  dz 

c 


j  (IV)  cos  - —  < 

_  4AmTim<»c(~I)'"'*''*(l  +<<.i)tanh-rin,c 

'  ™  (yL  +  (nit/cp)’ 


J  -J 

cos - dx  dz 


=E 


^^nnacl-l)"* '*•'(14- ^onltanh a„na 
Omni’  (aj,„  -(-  (I’T/o)') 


4- 


(ne6/c)^(Iff/o)^ 
Omn  4-  {li’la? 


rf 

4  e  J  —  a 


(VI)  cos  — —  cos  dxd: 

a  c 


_  B„lbc{  1  +  Sol){l  +  ion) 


2W’»»4-  f- 


-rio 


I3„,6cothl3„|6)  4-  I3„i6tanh0„|6] 


^oi  = 


on  — 


{. 

{ 


when  I  =  0 
when  /  7<^  0 

when  7?  =  0 
when  n  /  0 


where 
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The  individual  integral  components  may  now  be  summed  to  give  the  following. 


4A„„a^„ 6^)  tanh 


(«J.„  +  (Wa)’)" 

4An,.iac(-l)"‘'*'‘(l  +  ^tmltanhomaa  j  j  {nrb/c)^{lr/ay 
(ai.„  +  (Ir/ap)  v  +  ^  {Ir/a}^ 


B„iab(l  4-  ^of)(l  +  ^<m) 


Bnlbc(  1  +  Sot)(l  +  ) 


(lire \  ^ 

-^j  [(1  -  |0„|icoth/3„j6)  + /3„ifctanh/3„|l)] 
+  (~~)  ((1  -  /3„i6coth /9„|fc)  +  /3„|fctanh/3„i6] 


4D|„ttc(-l)"‘'*''*(l  +  <oiHanfaliimC 

-rimi’  (7L  +  {Wc)’) 


.^»,2_2  4.  (ivk/ul’ln'/c)’ 

i/m  IT  +  - 


.  4D|„7f„ac(-l)"'+"(l  +  «^)t»nh7,„c 
'  (Tj^n  +  (nfr/c)*)' 


Simplifying  this  equation  gives  the  following  result,  as  presented  in  equation  (13). 


4A,„„(-l)'"'*-*(l  +  <„„)tanhan,n° 
l’{al,„  +  {lir/af) 

4Am(-l)’"^"(l+<><,i)tanh7i„.c 

>’{fL  +  (n^/';)^) 

+  az"' 7  +  ^<m)('  +  +  c’l’) 


=  E- 


+  L' 


X  (2i/+  1  -/9„|6cothA,i6  +  /3„i6tanh/3„|6) 
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